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Abstract 

We investigate the possibility of a bosonic preheating in the simplest model of 
super symmetric Hybridinfiation (F-term inflation), which was considered first by 
Dvali et al. Here the inflationary superpotential is of the O'Raifertaigh-Witten 
type. The end of inflation is related to a non-thermal phase transition, which in 
the context of left-right symmetric models lowers the rank of the gauge group. 
Using the homogeneous classical field ansatz for the appearing condensates, our 
results indicate that the parametric creation of bosonic particles does not occure 
in the model under consideration. 
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1 Introduction 



The recent success in detecting neutrino masses has caused a renaissance of 
the idea of grand unification. In left-right-symmetric models, such as SU{3)c x 
SU{2)l X SU{2)r X U{1)(^B-L) or 5*0(10), massive Majorana neutrinos arise nat- 
urally. As the difference of baryon- and lepton number is a gauged symmetry 
in these models, there could exist interesting mechanisms for the creation of the 
corresponding asymmetries. On the other hand, supersymmetry not only pro- 
vides the right mass scale for the heavy neutrinos, but also gives the possibility 
of a inflationary potential in the GUT-Higgs-Sektor. Then it is possible to con- 
struct a cosmological model based on a consistent supersymmetric grand unified 
theory, which can be judged by both, its cosmological and its particle theoretical 
features. In this paper we will concentrate on an O'Raifertaigh-Witten model, 
which by the autors of was shown to be a realisation of Linde's Hybrid In- 
flation scenario p|. The superpotential can serve as a part of the Higgs sector 
of a left-right-symmetric model with global or local supersymmetry. The phase 
of inflation and the formation of density fluctuation in this model have been 
studied in 0]. Since a successful cosmological model needs a very effective 
mechanism for particle production after inflation, here we investigate the possi- 
bility of a bosonic preheating in the inflationary potential under consideration. 
As preheating possibly allows for the creation of superheavy particles, it is a 
very interesting scenario not only in the context of lepto- and baryogenesis but 
also for the creation of dark matter. Focusing on qualitative insights rather than 
quantitative accuracy, we tried to work with a minimum of numerical expense. 
Our investigation was inspired by those in |^ , where preheating in the original 
version of Hybrid Inflation was studied. Our numerical results, which originally 
were presented in the context of an investigation of inflationary supersymmetric 
S'O(10)-models in P], are very similiar to those in a recent paper, based on a 
version of a NMSSM But, as shall be seen, we disagree with the physical 
interpretation of the results given there. In contrast to [0, we do not see any 
evidence for bosonic preheating in supersymmetric hybrid inflation. 

2 Supersymmetric Hybridinflation: the Basic 
Scenario 

A common way to achieve an inflationary scenario from a supersymmetric theory, 
is to use the simplest O'Raifertaigh-Witten model for spontaneous symmetry 
breakdown It is given by the superpotential 



where the couplings are protected by either a continouus or a discrete i?-symmetry, 
which could descend from a string theory. The superfield C = (ac'ipc) is a rep- 
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resentation of the Lie-algebra Ac of the Lie-group G, and C = {dQ,ijjc) is the 
conjugate representation. Here the a's represent the scalar component while 
the ijj^s are the spinors. The Superfield X = [ax.ipx) is an Ac-Singlet. The 
mass-scale /i could be caused by a string compactification. In the case, that G 
represents a gauge symmetry we have the tree-level scalar potential 

-| dim(G) 

V^'^ = \Fx\' + \Fc\' + \Fa\' + ^ E W (2.2) 

^ r=l 

dim(G) 

= K^laxacl"^ + K^laxacl"^ + \K,acac - fi"^]"^ + - X! l^'A'^- 

Here represents the vacuum expectation value (vev) of the quantum field 
(li. Then V^^^ is minimized for argag + arg ac = and it is independent of 
argac + arg ax and arg ac + arg ax- Supersymmetry requires the D-terms to 
vanish identically for each single generator of the group G. If G represents a 
product group including an f/(l)-factor, this means that lad = |ac|. 

In the following we will be interested only in the direction of the Ag-multiplets 
ac and a^, which aquire a G breaking vev. Therefore they will be denoted by the 
same names. Now, concentrating on the D = 0-direction, there exists a certain 
/^-transformation which brings the scalar components a^ to the real axis. The 
corresponding canonically normalized scalarfields 0, a and couplings A, g are 
given by g]: 



or := 2ac = 2ac (2.3) 
(f) := V2ax (2.4) 

A := ^ (2.5) 

, := ^ (2.6) 
M ■= (2.7) 

where typically = 2A. Using this conventions the scalar effective potential in 
terms of the fields and a reads 

V{^, a) = ^ (m^ - Aaf + + V«(0), (2.8) 

which resembles very much Linde's original model of Hybrid Inflation 0. Here 
V^^^cj), a) represents the loop corrections to the tree level potential, which vanish 
for the supersymmetric case. For a=0, (p > M/g =: (pc the gauge symmetry 
remains conserved while supersymmetry is broken in the sector of the gauge 
singlet X inducing loop-corrections V^^\(l),a). Including 1-loop-corrections the 
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effective potential in this regime is given by 



V{<P > 0,) 




where A is a renormahzation scale. Since soft susy breaking terms lead to scalar 
masses of 0{TeV), such terms are neglegible compared to the GUT-scale mass 
parameter M. For a large range of the parameters k and fi the potential above 
satisfies the Slow Roll conditions for inflation. Here inflation is caused by the 
"cosmological constant" The scalar gauge singlet is the only degree of free- 
dom, which then has a nonvanishing vev 0. This vev very soon is dominated by 
the zero momentum Fourier mode and it forms a condensate with a very small 
extension in momentum space. Referring to [y] it can be approximated by a ho- 
mogeneous classical scalar field. For this reason we will call the inflaton and 
identify it with the zero momentum mode. The slow roll regime of the singlet is 
a typical aspect of the superpotential (|2.1|) and has become popular as effect of 
"a sliding field" outside the context of inflation [|1^ . 



The end of inflation is connected to the end of the gauge symmetric phase: 
as soon as reaches the critical point 0c the effective mass = — + (7^0^ of 
the Higgs-field vanishes. The classical equations of motions alone cannot tell us 
something about a phase transition, since the derivative \^(0c, a = O),^- vanishes 
identically and the vev of a remains zero. But, following the argumentation 
of Garcia-Bellido and Linde [Q, as slides towards zero, quantum fluctuations 
around this vev will get tachyonic masses stimulating an exponential growth of 
modes, whose momenta are smaller than the effective mass, k < |mo-|, where k 
means the comoving momentum. The modes with k > |mo-| will not grow at 
all. The result is a spontaneous breakdown of the gauge symmetry caused by the 
inhomogenous distribution of the field a with (a) = 0. Inflation ends typically 
near this phase transition. 

We now want to understand, if the behaviour of the Higgs field a can be 
approximated by an homogeneous classical field. As pointed out in [^] the dis- 
tribution of the Higgs field would be homogeneous on scales / ~ or even 
somewhat greater, if was constant in time. However, we will see soon that 
this mass is oscillatory. This changes the range and the mechanism of ampli- 
fication of the "tachyonic modes". Due to this fact, the distribution crfc(t) in 
momentum space will be even more localized around the zero momentum mode, 
such that the spatial distribution will be homogeneous even on a scale / which 
is somewhat greater than Imo-I^"*^. At least at such scales a{t,x) also can be ef- 
fectively described by a classical homogenous scalar field cr(t), i.e. as the zero 
momentum mode of such a classical field, which rolls down from the critical point 
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to the supersymmetric minimum, where a = M/\/\. We will call a{t) the Higgs- 
condensate. The range of validity for this ansatz will become clearer during our 
numerical investigation. 

As pointed out in local supersymmetry modifies this scenario because 
of the non-renormalizable terms, whose influence on the inflation depends on 
the value of the coupling constant k. But in most of these cases the end of 
inflation is dominated by the renormalizable terms in the potential. Thus, non- 
renormalizable couplings will not contribute after the phase transition and will 
be ignored in the following investigation. 



3 Evolution of the Background 

After the phase transition both of the fields x) and x) form condensates, 
which essentially can be described by the two homogeneous classical fields (pit) 
and cr(t). Free Higgs or singlet particles will be identified with quantum fluctua- 
tions around these condensates. In this section we concentrate on the evolution 
of the background and we neglect the effect of the quantum fluctuations. For 
reasons of simplification we use a "M-rescaling" to natural variables: 

t y:= Mt 
X C, '■= Mx 

k ^ K:=-^ (3.1) 



a ^ s :- 



M 
M 



M Or. 



Here the comoving momenta K and the scale factor a(?/) are normalized in such a 
way that ci(?/c) = 1 is fullfilled, where means the M-rescaled time at the phase 
transition. 

Neglecting the sub-dominant loop corrections, which certainly will appear, as 
long as the supersymmetric minimum is not reached, the equations of motion 
after the phase transition read: 

(^ + 3%)4 + 23^(y))/(y) = (3.2) 
+ + (-1 + f\v) + s\y))) s{y) = (3.3) 

where h is the M-rescaled Hubble Parameter, h := {^a)/a, with 

h\y) = [2f (y) + s'iy) + (1 - s'{y)r + 2f{y)s\y)] . (3.4) 

pi 
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Figure 3.1: Evolution of the back- 
ground 
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Figure 3.2: Oscillations in the 0-cr- 
plane 



Important qualitative insights then are possible without any detailed calculation: 
Both of the classical fields, / = ^ and s = effectively oscillate with the 

same frequency, approximately given by = -\/2, around the supersymmetric 
mmimum, = 0, a = M/yX. For energetic reasons the phase difference amounts 
to 7r/2. Coming from the critical point (0 = M/g,a = 0), depending on the 
initial values the dynamics in the — cr-plane should tend to be nearly one 
dimensional. Since in our case the initial time derivates of the fields should be 
very small slow roll values the trajectory should be near to a straight line. Another 
interesting point, obvious without any calculation, is that the dynamics of the 
post-inflationary system, measured in its natural time and lengthscale M~^, only 
depends on the scale of the phase transition, M/ ^/X. As the M-rescaled classical 
fields / and s and their derivatives, take maximal values of 0{1) effectively, this 
scale has influence on the damping but not on any other aspects of the background 
oscillations: the smaller the scale is, the longer the system swings correponding 
to its natural timescale M~^. Thus, qualitative results of one breaking scale will 
also appear at another breaking scale. 

This discussion is confirmed by our numerical integration of the system of 
equations, as shown in figures and |3]^. The parameters used in this and all 
of the following investigations are given by 



A = 0.625 X 10"^ 

/ = 0.125 X 10"^ 

M = 0.350 X 10^^ (3.5) 
M 1. 

= 1.40 X lo^^ 

which corresponds to a phase transition near to the SUSY-GUT scale in a super- 
symmetric 50(10) model in 0. The solutions are essentially damped oscillations 
around the sypersymmetric minimum, with the predicted phase. The dynamics 
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Figure 3.3: energy densities of the background 



R 


■■= pa^ 


R<j, 


:= p^a^ 


Ra 


■■= Pad^ 



in the /-s-plane is practically one dimensonal and can be restricted to the straight 
line for very small initial values of the time derivatives of the fields. Anharmonic 
behaviour is caused by the interaction of the classical fields and by the fact, that 
the Higgs field has a negative mass contribution from the tachyonic mass. 

Our result for the supersymmetric case {g"^ = 2A) is quite different to the case 
with (yf^ = A , which was considered in 0], where the classical fields oscillate in a 
rather chaotic way. But it essentially equals the result in 

Now we consider the comoving energy densities after the phase transition: 

(3.6) 
(3.7) 
(3.8) 

Here approximately p = M^/AX is the system's total energy density at the phase 
transition. As the Universe passes a phase of matter domination after inflation, 
such that p ~ a~^, the product pa^ will be constant. This is what our numerical 
investigations show in figure |3.3| . After the phase transition R is growing for 
a short time, up to the value of 1.24M'^/4A and then remains practically con- 
stant. It takes the background about 120 oscillations until it swings in such a 
,, harmonic" way, that the universe behaves matter dominated. Since we were 
interested in the effective distribution of energy between both of the swinging 
classical fields, in the calculation of R^ and R^ we integrated out the short time 
dynamics {t ~ 0{M~^)) by averaging over the the background oscillations. Here 
we respected the fact, that the frequency of the oscillations is slightly varying with 
time. The plot shows, how energy is exchanged between the classical fields - one 
reason for anharmonic behaviour in figure After about 2000 oscillations the 
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energy distribution between the oscillation of the inflaton and Higgs-condensate 
is approximately 1/3 : 2/3. Then the classical fields essentially are decoupled 
oscillators. 



4 Parametric Excitation of Quantum Fluctua- 
tions 

Our investigation of the quantum fiuctuations uses the M-rescaling introduced in 
section ^. We studied the mode equations following from the Heisenberg expan- 
sion of a quantum field which stands for the M-rescaled quantum fiuctuations 
6s{y,^) and Sf{y,^) around the classical fields s{y) and f{y)- 



(27r)3G3/2(^) 



(4.1) 



where p = K/a is the M-rescaled physical momentum corresponding to the M- 
rescaled commoving momentum K. The operators fe^'*'^ satisfy the commutation 
relation: 



Using the co moving momenta K = \K\ = |p|a for notation, this reads 



1 



{2nf 6{K - K'). 



(4.2) 



(4.3) 



Then the occupation number is given by 



riK 



UJ 



+ \Xk\' 



K 



and the particle density reads: 



(27ra)= 



tHk- 



(4.4) 



(4.5) 



Neglecting higher order terms in the quantum fiuctuations, the mode equations 



dy' 



+ 



+ 2s' - 



SMv) = 



+ i — + 3s^ + f-l- 5sKiy) = 0, 



dy 



(4.6) 
(4.7) 




As numerical calculations show, very soon (3/2)(aa + a?' /2)/{a?) ~ is fullfilled 
within the numerical accuracy and this term can be neglected. This corresponds 
to a matter dominated universe, where typically aa = — . We performed a 
numerical integration of the full set of differential equations. The values of K, for 
which these investigations were done, were selected by a pre-investigation using 
the well known approximation by the Mathieu equation |Tl[ . 



Before turning over to our numerical results, it is sensible first to consider 
the problem of particle interpretation. As clearly pointed out before, both of 
the classical fields, and a, are approximate representations of the quantum 
fields 4>{t,x) and a{t,x) in the low momentum zone. Although we identified 
these classical fields with their zero-momentum modes in principle they have an 
extension in momentum space. The edge of the condensates in momentum space 
will be situated at a momentum which is somewhat smaller than the mass of 
the corresponding field. In this sense the use of a homogeneous classical field 
is similiar to a limitation of the ,, resolution" in an optical investigation. But if 
this is true, interpreting the quantum fluctuations around these classical fields 
as free particles in presence of a homogeneous background is justified only for 
those modes, whose momenta are well above the edges of the condensates. This 
important point was apparently not taken into consideration by the authors of 

0- 

As in the last section, we performed numerical calculations of the full system 
of differential equations mentioned above. We did not find any parametric am- 
plification of the Higgs-fluct nations with K > 1/3. For i^' < 1/3 there exists an 
amplification which can be described as a mixing of both effects: mainly para- 
metric resonance and partly tachyonic mass. In the well known picture of the 
Mathieu equation, 

il? + ^fc - 2g cos(2z) j Xk = 0, (4.8) 
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which is apphcable restricting only to some few oscillations, this resonance starts 
inside the first resonance band {1 — q < Aj^ < 1 + q) with a resonance parameter 
q ~ 0.9, which corresponds to a strong narrow resonance regime. The amplifi- 
cation of these modes is essentially independent of the momentum and is given 
by the solution for K = 0. Furthermore, these modes lie in a range of mo- 
menta, which we estimated in section to have non-vanishing fourier modes of 
the Higgs-condensate: Because of this we do not consider these soft modes to be 
free particles. We think, that using the homogeneous classical field ansatz it does 
not make sense to study Higgs fluctuations with such momenta. Our result only 
tells us someting about the extension of the distribution of a which was described 
in section 0. 

The picture in the case of inflaton fluctuations is similiar. Here a tachyonic 
mass is missing and the effect of the resonance is much weaker. But there is no 
amplification of modes with K > 1/10 and similiar to the above case,the observed 
weak parametric amplification for modes with K < 1/10 cannot be intepreted as 
a creation of free particles. 

In order to learn more about the space/momentum dependence of the con- 
densates one could calculate the mode equation of the fourier modes sk and fx 
of the condensates. Since the original equations of motions are non-linear, the 
mode equations then are be rather complicated integro-differential equations. A 
first step of simplification would be to let their masses only be given by the zero 
momentum modes of the condensates. As can be checked easily, the mode equa- 
tions for fx are exactly the same as for dfx- Those for sk will look a little bit 
different from those for Ssk but will have similiar resonance parameters in the 
Mathieu approximation. Therefore, from the investigation that we did for the 
quantum fluctuations one can conclude that there will be a similiar resonance 
effect in a similiar range of momenta for the modes sk and fx- 

This shows that after the phase transition the inflaton condensate slightly 
looses its homogenity, while the Higgs condensate is non- homogeneous from the 
very beginning. However, at a physical scale which is somewhat bigger than 
aM~^ they are essentially homogeneous. Quantum fluctuations around these 
condensates should have monumenta which are bigger than the inverse scale of 
homogenity, in order to have a particle interpretation. 

Until now our investigations did not consider the effect of backreaction of the 
quantum fluctuations on the classical fields, which is very important in the case 



of a preheating with a very efficient particle production. As pointed out in |]lT 



the main effect of the backreaction is to lower the resonance parameters in the 
picture of the Mathieu equation. This means, that appearing resonances will 
become weaker and will last for a shorter time. Then, in our case, the general 
picture is not changed. The extension of the condensates in momentum space 
probably will be smaller, such that they will be even "more homogeneous" than 
without backreaction. But there will be no possibility of having a parametric 
amplification of fluctuations, which were not amplified also without backreaction. 
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Figure 5.1: parametric excitation of the k = M-mode of a scalar field x 

So we conclude that there does not exist a creation of free Higgs- or singlet- 
particles by parametric amplification. This corresponds to the result in for 
the different case of X = g'^. 

5 Parametric Excitations of External Fields 

Although we did not find a preheating of the quantum fluctuations, there could 
be a parametric creation of an external scalar quantum field x, that couples to 
the Higgs sector via 

+ IhWx'. (5.9) 

Then the equations of motions will be modified in an obvious way. The mode 
equation for the modes XK{t) reads 

( (P ( ^ ,2 hi ,\\ , , , 

(^+(^ + ^/^ + f.^]jx.(.,e)-0, (5.10) 

where again da = — was assumed. 

We used the same methods as before in our numerical investigation. Our 
result is very similiar to the result in for g'^ = A. We found a strong resonance 
only for the case without coupling to the Higgs-condensate, h2 — (see figure 
|5.1| ). But this situation is impossible, if we use the superpotential ( p.l[ ). The F- 
terms always lead to /12 7^ if /ii 7^ and vice versa. As can be shown very easily 
[Q, in this model there does not exist any coupling allowed by supersymmetry 
and leading to a bosonic parametric resonance. This result, as far as we can see, 
buries the hope for the existence of a bosonic preheating. 
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6 Consequences for Left-Right Symmetric Mod- 
els 



In left-right-symmetric models the phase transition that leads to the gauge sym- 
metry of the Standard Model, necessarily lowers the rank of the Lie group. The 
only renormalizable superpotential with this feature has the structure 

W = XCC + polynomial in X, (6.1) 

where X is a gauge singlet and C, C are appropriate spinor representations of 
the left-right-symmetric gauge group. It is an interesting question, if, using the 
superpotential (p.l|), there could be a natural embedding of inflation into the 
context of a right-left-symmetric model such as S'O(IO) or S'f/(3)c x SU{2)l x 
SU{2)ji X U{1)(B-L), which leads to a senseful cosmological model. 



In minimal S'O(10)-models [0, |T3| for example C should be a 16- dimensional 
spinor representation, the scalar component of which aquires a GUT-scale vev 
in the S'f/(5)-singlet direction during the phase transition. Since the SO (10) 
(more precisely: Spin{10)) is a simple and simply connected Lie group, breaking 
it down to the SM-gauge group will lead to unwanted monopoles. They appear 
in a first phase transition, when a 45 dimensional Higgs representation aquires a 
suitable vev to break Spin{10) down to the left-right-symmetric group Glr '■= 
SU{3)c X SU{2)l X SU{2)r x U{1)^b-l)- The inflaton then dilutes the unwanted 
remnants and ends in a phase transition which breakes the remaining gauge 
symmetry to the standard model group. This picture remains the same in a pure 
Giij-Model, C being a SU{2)r x t/(l)(B-L)-Higgs-doublet. 

Our investigations indicate, that in such models free bosonic particles cannot 
be created by the effect of parametric resonance after inflation. The parametric 
amplification in the ,,low momentum zone" of the quantum fluctuations, which 
appears in our numerical investigation, cannot be interpreted as a production of 
free particles. In this area not only the dynamics of the mode functions are es- 
sentially independent of the momentum but also the validity of the homogeneous 
classical field ansatz is doubtful. In our view, these numerical results perhaps can 
tell us something about the extension of the condensates in momentum space but 
shurely show us the limitation of the homogeneous classical field ansatz. 

Supposing that bosonic preheating does not take place, a next step would be 
to study the possibility of fermionic preheating for this superpotential, which re- 
cently was shown to be possibly very efficient . Since we deal with a supersym- 
metric theory, there is no reason for neglecting the fermionic superpartners. Also 
a parametric creation of Majorana neutrinos could appear, which should be very 
interesting in the context of leptogenesis. In this context it also would be nesse- 
cary to rule out the creation of the helicity-l/2-gravitinos by non-perturbative 
effects, which recently was shown to be of possible danger for cosmological models 



which involve supersymmetry []T3[. If fermionic preheating turned out not to be 
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efficient, particle production in left-right-symmetric models could be caused only 
by perturbative effects. Following [|l^ then still the production of super heavy 
matter could be possible. But, as we work with two condensates, both of them 
would have to decay very efficiently, in order to be not in conflict with standard 
cosmology. Considering the gravitino constraint on the reheating temperature of 
the Universe, this would mean strong limitations to possible couplings to other 
fields. Concluding, there is still some work to do, until the particle creation after 
inflation in this very simple supersymmetric model will be fully investigated. 
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